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The colloidal gel and glass transitions are investigated using the idealized mode coupling theory 
(MCT) for model systems characterized by short-range attractive interactions. Results are pre- 
^ l' sented for the adhesive hard sphere and hard core attractive Yukawa systems. According to MCT, 

, the former system shows a critical glass transition concentration that increases significantly with 

' introduction of a weak attraction. For the latter attractive Yukawa system, MCT predicts low tem- 

perature nonergodic states that extend to the critical and subcritical region. Several features of 
the MCT nonergodicity transition in this system agree qualitatively with experimental observations 
on the colloidal gel transition, suggesting that the gel transition is caused by a low temperature 
y j — ^ , extension of the glass transition. The range of the attraction is shown to govern the way the glass 

Q ■ transition line traverses the phase diagram relative to the critical point, analogous to findings for 

C/3 ' the fiuid-solid freezing transition. 

^ . 

, PACS numbers: 64.70.Pf, 64.75.+g, 82.60.Lf, 82.70.Dd 

a 

Ch ■ I. INTRODUCTION 

o 
o 

Attractions among stable colloidal particles lead to a diverse phase behavior. Colloidal attractions, unlike the 
molecular attractions, act usually over a relatively short (compared to the particle size) range. It is by now well 
established that when the range of the colloidal attraction is decreased the phase diagram undergoes a progression 
from gas-liquid-solid to fluid-solid coexistence, the latter with a subcooled critical point which is metastable with 
^vq respect to fluid-solid coexistence [Q-^ • Numerous experimental studies show that suspensions often form incompletely 
p\| ', equilibrated solids with the appearance of gels where one expects a fluid-solid |^~|ljl or gas-liquid ||l^ phase separation 
from equilibrium theory. The systems studied include mixtures of colloids and non-adsorbing polymer P pTl and 
sterically stabilized colloids in marginal solvents In the former case, the attractions stem from depletion of 

. the polymer coils from the regions_between closely spaced particles 18[ 1^] , and in the latter they are caused by surface 
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grafted chain-chain interactions [gO[ . The gel transition is observed when the range of attraction is short compared to 
the particle size. In the colloid-polymer mixtures this is achieved by choosing a small ratio of polymer to colloid size, 
whereas the overlap length of the surface grafted chains sets essentially the range of the attraction in the sterically 
^ , stabilized particle systems. 

■ While the equilibrium phase behavior of these systems is well understood the nature of the gel transition 

^ remains to be clarified. The gel state appears to be related to a ramified structure with interconnected particle clusters 
O [^,0. Temporal density fluctuations are very slow close to the gel transition and the suspensions acquire a yield 
^ stress and a finite low-frequency elastic shear modulus in the gel j|ll| , |l3| , p^ 17 2l|]. In the past the transition to the 

^^r|i|,|i|,p 



gel state has most often been interpreted as either a static percolation transition ]13| , p^ , p^ , where a sample-spanning 
cluster of particles forms, or due to the fluid-solid phase transition [Q. Comparison between integral equation 
^ , predictions for the percolation transition and experimental data, however, shows that the gel transition is confined to 
' the region in the phase diagram between the static percolation threshold and the gas-liquid spinodal ||lj,|l5| . Colloidal 
gels have also been attributed to states inside a gas-liquid binodal which is metastable with respect to fluid-solid 
coexistence. Such metastable binodals have indeed been observed for suspensions of globular proteins which also 
form gels when the ionic strength is sufficiently high ||23|-p5[|. 

In this work we present an alternative interpretation of the dynamical arrest of the gel structure which causes 
colloidal systems to become disordered solids. We propose that colloidal gels are nonergodic systems that form when 
a dynamic gel transition is traversed. We further suggest that the gel transition is a low temperature extension of the 
liquid-glass transition. The gels, however, differ physically from colloidal hard sphere glasses in that they generally 
display a larger elastic shear modulus and that the particles are more strongly localized in the gels; both of these 
effects are due to particle clustering induced by a short-range attraction among particles. 
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We demonstrate that this is a possible explanation for colloidal gel formation by applying the idealized mode 
coupling theory (MCT) of the form used to study the liquid-glass transition l2^-|3C|] to systems in which the attraction 
is restricted to short ranges. This study is motivated by the observations made by Verduin and Dhont [|l5|, who noted 
a structural arrest (nonergodicity) in connection with the gel transition similar to that observed for hard sphere 
colloidal glasses [pT|-p4| . Also Poon et al. have made such observations in the so-called transient gelation region 
of the colloid-polymer phase diagram for short polymers. Krall and Weitz |3^] have shown that, in the limit of 
strong particle aggregation, suspensions become nonergodic even at low colloid densities. To date, however, only a 
speculative connection has been made between the gel and liquid-glass transitions |p^ . 

We conduct a study of ergodicity breaking in two model systems: Baxter's adhesive hard sphere (AHS) system and 
the hard core attractive Yukawa (HCAY) system. Both systems supply analytical solutions for the static structure 
factor, the former within Percus-Yevick (PY) theory p6[ |, and the latter within the mean spherical approximation 
(MSA) [^-pot. This study provides more information on the AHS phase diagram and, in addition, serves to comple- 
ment a recent independent MCT study |Q on the temperature dependence of the AHS glass transition. Further, the 
HCAY system provides a likely candidate for the gel transition in colloidal systems as an ergodicity breaking dynamic 
transition of the same type as the liquid-glass transition, suggesting that the experimentally observed gel transition 
is a low temperature extension of the glass transition. 

In what follows, the idealized mode coupling theory of the liquid-glass transition, suitable for colloidal suspensions, 
is briefly summarized. Results for the temperature dependence of the glass transition are then presented and compared 
to the AHS phase diagram as predicted by density functional theory. Results are also shown for the HCAY system, 
which show that the MCT glass transition extends to the critical and subcritical region at low temperatures. The 
way in which the glass transition line traverses this part of the phase diagram is shown to be a strong function of the 
range of attraction. 



II. MODE COUPLING THEORY 



The mode coupling theory (MCT) of the liquid-glass transition provides a dynamic description of the transition 
p6|- ^ , ^ . For sufficiently strong interactions the dynamical scattering functions do not decay to zero with time, 
leaving instead finite residues - the nonergodicity parameters, also known as Edwards-Anderson parameters or glass 
form factors. Generically, concurrent with this long-time diffusion ceases and the zero-shear viscosity diverges, both 
due to a diverging relaxation time. This structural relaxation time is in turn related to the particles' inability of 
escaping their nearest neighbor cages. The glass transition within the framework of the idealized MCT is not a 
conventional thermodynamic phase transition; the constrained motion of the particles leads to a difference between 
time and ensemble averages, i.e. an ergodicity breaking transition. 

Application of the MCT to the Mori-Zwanzig reduced equations of motion for the density correlators, together 
with at^oo limit, leads to the following set of closed equations p6|-po|,p2[ 



^^(q, k) = q • (q - k) C|q_k| + q • k Cfe + qpc^^'> (q, q - k) 



(1) 



(2) 



where p is the number density, Sq is the static structure factor, Cg — {Sq — l)/pSq is the Fourier-transformed Ornstein- 
Zernike direct correlation function, and c^"^-* is the triplet direct correlation function. In this study we use primarily 
the so-called convolution approximation (c*-'^-' = 0) Jp7| , p9| j4^ for the triplet direct correlation function. Note that 
the coupling vertices ^(q, k) and V^''(q, k) in Eqs. |l|and |^ are free of singularities and vary smoothly with a set of 
external control parameters, e.g., the particle density. The coherent (/g) and incoherent (/g) nonergodicity parameters 
are defined as the long-time limits of the intermediate scattering function Fq (t) and the self-intermediate scattering 
function Fq (t) , according to 

fq= l\in{Fq{t)/Sq) (3) 

t — >oo 

/I-/™ (^.^W)- (4) 

The nonergodicity parameters determine the properties of the glass. The zero-frequency elastic shear modulus of 
the colloidal glass (in units of ksT/a^, with ksT the temperature and a the particle diameter) is given by p7| , ^ 
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The incoherent nonergodicity parameter is found to be well approximated by a Gaussian, the half-width of which 
is proportional to the mean-square displacement in the glass state p^ , p7| . The localization leiigth is defined as the 
root-mean-square displacement in the glass, and is determined from — I ~ q^rl for q ^ 

Eqs. |l| and ^ are solved self-consistently for the nonergodicity parameters as functions of specified external control 
parameters: the reduced temperature r and volume fraction (f> = irpa^ /& for the AHS system, and the reduced 
temperature K~^, screening parameter b, and volume fraction for the HCAY system (cf. below). The solution proceeds 
by iteration, first on /q, and subsequently on . Transition lines delineating ergodic and nonergodic states were found 
by bracketing, and the monotonicity property of the iteration |^ was employed. The integrations were performed 
numerically using Simpson's rule on a uniformly discretized wavevector grid: qi = i/S.q, i = Q, . . . , N . The parameters 
Ag and N used varied somewhat, but most results were obtained using 0.15 < Aqa < 0.30 and 600 < N < 1000. The 
iterative solution scheme was complemented occasionally by an algorithm to speed up convergence by using stored 
previous iterates |Q. We have also directly integrated the equations of motion ^] (adjusted to obey Smoluchowski 
dynamics). This yields the entire time-evolution of the density correlators Fq{t) and F^{t), the long-time limits of 
which were found to be identical to the solutions of Eqs. |l] and ||. 

The result = /| = is always a solution to Eqs. ^ and |^, implying that correlations among density fluctuations 
vanish for long times. At low densities this is the only solution, hence the system is in a fluid, possibly metastable 
fluid state. Above a critical volume fraction also non-zero solutions appear, which correspond to nonergodic glass 
states. The physical solutions to Eqs. |^ and |2| corresponding to the long-time limits defined by Eqs. || and ^, were 
identified by choosing the largest solutions for fq and P^jEII • 



III. ADHESIVE HARD SPHERE SYSTEM 




The AHS pair potential consists of an infinitely deep and narrow well located at particle contact, given explicitly 

by [H 

u{r)/kBT^ lim I In iM|z£) a<r<d (6) 

where r is a reduced temperature and r is the interparticle separation distance. The r — > cxd limit of the PY~AHS 
system corresponds to the PY hard sphere system. Using this as the starting point, the earlier MOT result for the 
hard sphere glass transition volume fraction 0c = 0.516 was reproduced. With more accurate hard sphere static 
inputs one obtains (pc = 0.525 p7| , ^ , ^ . Experiments locate the colloidal hard sphere glass transition at 0c ~ 0.58 
pl|-^, showing that the idealized MCT prediction for the hard sphere 0c is too low; this is presumably caused by 
the strong restriction of the modes of structural relaxation imposed in the MCT. 

The locus of critical glass transition points is shown in Fig. [l] as a function of the particle volume fraction and 
the reduced temperature r, with the shaded region denoting nonergodic glass states. Upon decreasing r the glass 
transition point moves along the line Bl in Fig. |^ to higher density, contrary to the findings for particles interacting 
via a molecular interaction potential of Lennard- Jones form |^ . We see that starting with a hard sphere glass and 
introducing a short-range attraction leads to an ergodicity restoring transition, provided the density of the isochore 
is not too high. The glass transition for r = 2 occurs at a volume fraction of 0.5527, significantly higher than the 
MCT hard sphere result. The strength of attraction needed to shift the glass transition to higher density is weak; 
a second virial coefficient mapping reveals that a reduced temperature r = 2 corresponds roughly to a O.OScr wide 
square- well with a depth of ^0.6 ksT. An examination of the static structure factors along the critical glass transition 
boundary Bl in Fig. ^ shows that they are not markedly different from those of hard sphere suspensions; however, 
subtle changes in the structure lead to significant changes in the critical glass transition density. 

At high temperatures we observe that the MCT predictions for the localization length follow roughly a Lindemann 
criterion, given by the hard sphere result « 0.074ct, in agreement with the results of the Lennard- Jones study 
Below r ~ 5, however, the Lindemann criterion is violated, the particles being now more strongly localized in the 
glass state. We note also that inclusion of improved triplet correlations in the manner of Barrat et al. ji^, who 
used an approximation due to Denton and Ashcroft for c^'^) [ pO[ , leads to a qualitatively similar phase diagram with 
boundaries shifted slightly to lower densities and temperatures relative to those shown in Fig. |l|. 

The AHS phase behavior has been the subject of several studies, most of them using density functional theory (DFT) 
[^-p5[ . Selecting the most recent one by Marr and Cast |5^ for comparison, who used the modified weighted-density 
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approximation (MWDA) ||5^, we find that the Bl glass transition is confined to the metastable region between the 
fluid-solid coexistence lines (see Fig. |^). One striking feature is that the Bl glass transition line from MCT and the 
MWDA freezing transition line track each other, the quantity A(f> — (f>c — (l>f, with (j)f the volume fraction at freezing, 
being nearly temperature independent. 

This result has interesting consequences for the diffusion constants at the freezing densities . When sufficiently 

close to the glass transition the long-time self diffusion coefficient assumes its asymptotic behavior governed by the 
distance to the glass transition singularity. The normalized long-time self diffusion coefficient has been found to 
exhibit universality along the fluid-solid freezing transition [s^ . The comparison made here shows that this condition 
may be related to a universality of the proximity of the freezing transition to the glass transition. At least, it suggests 
a deeper connection between MCT for the liquid-glass transition and DFT, a topic that has been explored to some 
extent ]59| ]. 

In following the glass transition line from high to low temperatures (the line denoted by Bl in Fig. |l]), we flnd a line 
crossing similar to that studied within schematic (g-independent) models p7|, 47| , |60| . This region in the phase diagram 
has been studied in detail recently by Fabbian et al. [Q. At the crossing between the Bl and B2 glass transition 
lines, it is B2 that determines the behavior of the physical solution as the noncrgodicity parameters associated with 
B2 are found to be always greater than those associated with Bl. Thus, along the B2 line bordering the fluid phase, 
fq for each q jumps discontinuously between and flnite values, and between smaller and larger finite values when 
the B2 line is traversed in the glass. The appearance of the B2 line is a result of an endpoint (cusp, A3) singularity 
|^,|0|j60[|6| in the AHS phase dia gram, where three solutions of Eq. Q for fq coalesce. This singularity appears as 
the termination point of the B2 transition line. It is connected to another bifurcation point with triply degenerate 
fq solutions located at lower temperature by the B3 transition line shown in Fig |^. Neither the low temperature 
endpoint, the piece of Bl between it and B2, nor the B3 transition line play a role in determining the glass dynamics, 
but show the connectivity among the bifurcation solutions of Eq. |^. 

The B2 line in Fig. |l| exhibits unusual properties. Varying the numerical parameters N and Aq, such that the 
maximum wavevector qmax = NAq changes, shifts the location of the B2 glass transition line and the high temperature 
endpoint in the phase diagram. Such a variation is not observed in connection with the Bl glass transition line. In 
addition, we were unable to identify a set of N and Aq such that the fq associated with the B2 glass transition 
decays to zero within the prescribed wavevector range. The results shown in Fig. |l] were obtained using N — 700 and 
Aqa = 0.2. 

It is possible that the anomalous behavior of the B2 glass solutions results from the atypical behavior of the AHS 
Sq in the large q limit caused by the singular nature of the AHS pair potential. The AHS Sq decays slowly for large 
q as Sq 1 + 20Ab sin {qa)/qa, where Ab is the solution of Baxter's quadratic equation We do not consider the 
behavior of the B2 glass solutions here further; instead, we show in the next section that the HCAY system exhibits 
glass transition lines that extend to low temperatures and densities. Several properties of these low temperature 
glasses are in qualitative agreement with experiments on colloidal gels. 



IV. HARD CORE ATTRACTIVE YUKAWA SYSTEM 



In this section we examine the effect on the glass transition of introducing a finite range of attraction via the HCAY 
system. The HCAY pair potential is given by 

{oo < r < a 

_J<_e-Hr/a-i) ^ ^ ^ (7) 

where the dimensionless parameter K regulates the depth of the attractive well and the reduced screening parameter 
b sets the range of the attraction. Using the MSA static structure factor |^^-^ as input, the MCT was solved for 
three different screening parameters: b = 7.5, 20, and 30. The progression of the glass transition can be traced from 
the (PY) hard sphere limit, corresponding to K=0, to lower temperatures in terms of the reduced temperature K~^. 

In Fig. ^ we show the gas-liquid spinodal curves, studied in detail by Cummings, Smith, and Stell [^,|2j. The 
critical temperature is sensitive to the range of the attraction, decreasing with increasing b. The spinodal curves are 
shown as indicators of where gas-liquid phase coexistence will occur, should there be a stable liquid phase. Included in 
the diagrams in Fig. ^ are the corresponding loci of glass transition points. At high temperatures and small values of b 
(b=7.5) they are relatively insensitive to the strength of the attraction, showing only a minor initial movement toward 
higher densities. Increasing the value of the screening parameter b (b=20 and 30), which decreases the range of the 
attraction, leads to a small initial increase of the glass transition density upon lowering the temperature, although this 
trend is not as pronounced as in the AHS system; subsequently, at lower temperatures the glass transition is induced 
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at increasingly lower densities. For b=7.5 and 20 the nonergodicity transition lines reach subcritical temperatures 
and approach the liquid side of the spinodal. For b=30 the nonergodicity transition line lies entirely within the fluid 
phase above the two phase region, and extends to subcritical temperatures at low densities. 

The MCT used here does not account for large concentration gradients and additional critical slowing of relaxations. 
As there is no small expansion parameter in MCT, it is difficult to ascertain when, upon approaching a critical point, 
this form of MCT should be replaced by a more complete theory, including a more sophisticated handling of the 
critical dynamics (see ref. | |6^ and references therein). At high temperatures wavevectors around the primary peak of 
the structure factor contribute the most to the mode coupling integrals in Eqs. |l|and|^. With decreasing temperature, 
on the one hand, the small wavevector structure in Sq leads to a stronger coupling on large length scales; on the other 
hand, the attractive interactions become of increasing importance on all length scales. The former effect, which can 
be expected to appear for all ranges of attractions, leads to nonergodicity transitions which trace the spinodal lines. 
These transitions will be discussed in the appendix as here the present MCT is least reliable because it does not 
include all relevant mode couplings and will not correctly describe the dynamics near the critical points | |63[ |. The 
latter effect, important for systems with short-range attractions, can be studied by an asymptotic analysis of the 
MCT equations and will be seen to dominate the low density glass transitions for large values of b. 

At low densities and in the limit of strong attractive interactions, the Ornstein-Zernike direct correlation function 
becomes independent of density. Specifying this to the MSA of the HCAY system, this limit corresponds to (p ^ 
and K — > oo. Considering the asymptotic limit 

0-^0 and K ^ oo , so that F — — g— — constant , (8) 

the MCT equations simplify because 5^ — > 1 follows. The nonergodicity transitions then occur at F = Fc(b), leading 
to the asymptotic prediction Kc oc 1/ 

For short-range attractions, in the limit of b ^ oo, a further simplification arises because the coupling constant F 
approaches a unique value at the transition, Fc 3.02 for b ^ oo, and the nonergodicity parameters now depend only 
on the rescaled wavevector q = qa/h: f^iq). The asymptotic transition lines are shown in Fig. ^ as the chain 

curves. We find excellent agreement with the MCT transition line for b=30 at low density, demonstrating that the 
low density nonergodicity transitions are not driven by the divergence of the small wavevector limit of Sq. Moreover, 
the asymptotic model is seen to capture the behavior of the full MCT transition lines - where present - qualitatively 
and even semi-quantitatively at higher densities. 

We further point out that Eq. ^ for the single particle dynamics and, thus, the incoherent form factors /| are not 
dominated by small wavevector variations in the static structure factor. Instead, fq and the localization length Vg are 
dominated by the small distance or large wavevector behavior of the liquid structure. In the asymptotic limit of Eq. 
^ this also holds for the collective particle dynamics and f'^{q) = is obtained, where both functions show rather 

large non-Gaussian corrections. 

For the system with b=20. Fig. ^ shows that the glass transition nearly meets the critical point (see also Fig. Q 
in the appendix). This aspect is in qualitative agreement with the behavior of the sterically stabilized suspensions 
studied by Verduin and Dhont jist . They observed a gel transition that traversed the phase diagram from high 
density and temperature to the critical point. The transition, which they refer to as a static percolation transition, 
was associated with a non-decaying intermediate scattering function and non-fluctuating dynamic light scattering 
speckle patterns. Thus, it has the expected properties of the ergodic-nonergodic dynamic transition predicted by 
the idealized MCT. Moreover, they were able to follow the transition into the unstable region inside the spinodal 
curve, showing that complete phase separation does not occur because of interference from the gel transition. We 
cannot extend the calculation of the MCT gel transition line into the unstable region because an appropriate Sq is 
not available and the theory assumes closeness to equilibrium ||29| ]. 

Restricting the range of the attraction sufficiently, as for b=30. Fig. ^ shows that the glass transition line passes 
above the critical point and reaches subcritical temperatures on the gas side of the metastable spinodal. For such 
systems we may speculate that the glass transition renders the entire spinodal curve and the liquid phase dynamically 



inaccessible. This feature appears to agree with some measurements on sterically stabilized suspensions 14 |r^, in 
which only a liquid-gel transition was observed. 

Recent measurements by Poon et al. |0 suggest that the colloid-polymer mixtures with a small polymer/colloid 
size ratio (^ « 0.08) may belong to the class of HCAY phase diagrams with b < 20, where the nonergodicity transition 
line meets the spinodal on the liquid side. This interpretation includes a possible explanation for the growth of the 
small-angle scattering peak for samples with low colloid concentrations , and that the denser colloid domains arrest 
in the transient gelation region. 

To further clarify the physical mechanism of the gel transition and the properties of the gel states, various aspects 
of the solutions of the MCT will be discussed for the three cases: h—7.5, 20, and 30. In Fig. O we show the evolution 
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of the coherent nonergodicity parameter fq along the critical glass transition boundary corresponding to b=30 in 
Fig. ||. As seen, the width of fq increases with decreasing temperature (increasing K). This behavior of fq with 
decreasing temperature is a result of a corresponding increase in the range of Sq, which results from particles being 
strongly correlated near contact, i.e. due to particle clustering. For longer range attractions, like the b=7.5 system, 
the width of fq and fq remain essentially unchanged along the glass transition boundary, which reflects the lower 
degree of particle clustering in this system. Note that fq becomes a density and temperature independent function, 
fq — > f(qa/h), in the limit of strong attractions. This prediction is shown in Fig. ^as the bold line, and agrees almost 
quantitatively with the full MCT fq solutions for large values of K. 

The localization length in the glass decreases along the glass transition boundary when the attraction is sufficiently 
short range. This decrease in the localization length is shown in Fig. |[ and is caused by the increased contribution 
from large wavevectors in the MCT integrals in Eqs. |^ and For longer range attractions, like the b=7.5 case, the 
localization length stays close to the value dictated by the Lindemann criterion and found at the glass transition in 
the hard sphere system 26 i^]. Thus, for short-range attractions the particles are more strongly localized in the glass 
than for systems with somewhat longer range attractions. At low temperatures the localization length saturates at a 
limiting value which is inversely proportional to b because of Eq. ^, which leads to the prediction Vg — s- O.Olfi/b for 
b ^ oo. 

In addition to an increased width at low temperatures, the small wavevector behavior of fq changes dramatically, 
such that at low temperatures the intermediate scattering function for small q practically does not decay with time at 
all (see Fig. ^). This indicates that large scale assemblies of particles behave essentially as static objects, where the 
single particles are tightly bound to the particle clusters (see Fig. Q). The asymptotic limit, Eq. || and b cx), which 
results in /(g) — 1 (x q^, stresses that this is caused by the short-range attraction. Such a rise in fq for small q is 
observed also in the b=7.5 system, where it is caused by a different mechanism, namely the increase in the isothermal 
compressibility on approaching the gas-liquid spinodal. There, this leads to coherent nonergodicity parameters which 
are essentially hard sphere-like except for a large q ^ value. 

In Fig. 1^ we show the zero-frequency elastic shear modulus as a function of the reduced temperature along the glass 
transition lines in Fig. When the range of attraction is comparatively large (b=7.5) the shear modulus remains 
constant at the hard sphere value, even for suspensions near the spinodal. This illustrates that the shear modulus, like 
the localization length, is determined by the large wavevector behavior of the liquid structure Sq, and is unaffected 
by long- wavelength density fluctuations in our calculations. For shorter range attractions (b=20 and 30), the shear 
modulus is dominated by particle clustering; it increases strongly with decreasing temperature because of the stronger 
binding among particles, eventually showing a maximum for suspensions close to the bend in the versus (j)c curves, 
where (jjc begins to decrease strongly with decreasing temperature. At lower density the shear modulus becomes linear 
in the density at the transition, according to G oc K^(/)c/b, as predicted by the asymptotic solution in Eq. ^, and 
as observed in the dilute limit of the b=30 system. These results show that low temperature nonergodic structures, 
proposed to be colloidal gels here, are distinct from colloidal glasses in that they generally display a larger static shear 
modulus and strongly localized particles bound in clusters. 

To more clearly connect this study of the low temperature behavior of the glass transition to the experimental 
studies of the gel transition, we have calculated the intermediate scattering function upon approaching the glass 
transition at fixed volume fraction. This mimics the Verduin and Dhont study |]l5| , in which they performed low-q 
dynamic light scattering experiments on a series of suspensions at flxed volume fraction close to the gel transition. 
As noted already, the HCAY system with b=20 exhibits a glass transition line that nearly meets the critical point. 
This qualitative aspect is shared with the experimental phase diagram determined by Verduin and Dhont. We have 
selected four suspensions with b=20 and (j) = 0.4 at different reduced temperatures (shown as open circles in Fig. ||), 
such that the suspension with the lowest temperature (K=:12) is located in the glass. 

The resulting normalized intermediate scattering functions corresponding to these suspensions are displayed in Fig. 
^ for a fixed wavevector qa=0.2, the same wavevector as that used in the measurements by Verduin and Dhont. 
As the normalized intermediate scattering function is the quantity that one measures in dynamic light scattering 
experiments, we can compare Fig. ^ with the results of Verduin and Dhont (see their Fig. 11). This comparison 
shows excellent qualitative agreement between their dynamic light scattering data and our calculated Fq{t)/Sq. Away 
from the transition the decay of Fq(t)/Sq is approximately exponential for this wavevector. The decay becomes 
slower when the temperature is decreased until K = 12, when Fq(t)/Sq no longer decays to zero. Instead, a long- 
time plateau with a value near unity is obtained, which corresponds to the nonergodicity parameter fq at qa = 0.2. 
Note that additional incoherently scattered light due to particle size polydispersity Q may contribute appreciably 
and cause fq to attain such a large value. Nevertheless, the dynamical arrest of Fq(t)/Sq agrees with our proposed 
ergodic-nonergodic transition for the gel transition and is captured by the idealized MCT. 

As has been shown in the past, the range of the colloidal attraction dictates where the fluid-solid freezing transition 
passes through the phase diagram and whether there is a stable liquid phase pl-pt]. In the same manner, Fig. 
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shows that the range of the attraction determines how the glass transition traverses the phase diagram relative to the 
critical point. We have not compared the diagrams in Fig. y with results for the fluid-solid and gas-liquid transitions 
(see e.g., P, |6^ , |65[ ) . The MCT relies on the static structure factor input, which was provided using the MSA. For 
short-range attractions the MSA produces relatively poor structural and thermodynamic predictions. Thus, a fair 
comparison should be made with a theory, e.g., MWDA js^, which can use the same input as that supplied to the 
MCT. Alternatively, the MCT can be solved using a more accurate static structure factor, such as that from HMSA 
theory . This would enable a determination of the location of the glass transition relative to that of the fluid-solid 
freezing transition, testing the conjecture made here that the glass transition line tracks the freezing transition at 
higher density in the phase diagram. 

V. DISCUSSION AND CONCLUSIONS 

The idealized MCT has been shown to provide a possible explanation for important aspects of the colloidal gel 
transition. In this scenario the arrest of the dynamics during the gel transition is caused by a low temperature liquid- 
glass transition. The underlying phenomenon is a breaking of ergodicity, caused by long-time structural arrest. It is 
accompanied by the cessation of hydrodynamic diffusion and the appearance of relatively large finite elastic moduli 
as the particles are tightly localized in ramified clusters. 

The AHS system was found to have endpoint singularities in the phase diagram. The spectacular dynamics close 
to the MCT endpoint singularity has been the focus of a recent study However, the glass transitions that 

occur at low temperatures in the AHS system are accompanied by numerical difficulties, resulting from the singular 
nature of the AHS pair interaction potential. Nevertheless, the AHS system provides insight into the temperature 
dependence of the glass transition in systems with weak short-range attractions. Subtle changes in the structure 
caused by the attractions lead to an initial increase in the glass transition density with decreasing temperature. The 
particles forming the glassy cage tend to stick together, thereby creating openings in the collective cage around a 
central particle which have to be filled by increasing the critical colloid density. We suggest that the recrystallization 
of glass samples at high densities upon addition of short polymers, reported in is explained by this shift of the 
glass transition density to higher values. Moreover, the MCT glass transition line was observed to lie parallel to the 
DFT fluid-solid freezing transition at high temperatures in the AHS phase diagram. 

Introduction of a finite range of attraction and replacement of the FY theory with the MSA via the HCAY system, 
yields glass transition lines that extend to low temperatures in the phase diagram. For HCAY systems with moderate- 
range attractions the glass transition line crosses the liquid binodal. When the range of the attraction is further 
restricted the glass transition line passes above the critical point, likely rendering part of the (metastable) equilibrium 
phase diagram irrelevant. Preliminary solutions of the dynamical MCT equations for the b = 30 HCAY system 
indicate that nearby Ai singularities with I > 2 appreciably distort the time dependent structural correlators in the 
intermediate time windows, even though no singularity could be found in the phase diagram. 

The nonergodicity transitions of the HCAY system are influenced by two mechanisms which are absent or not 
dominant in the hard sphere and Lennard- Jones p9| systems. In the latter two, the temperature dependence of the 
critical density pc is either trivially absent or arises from the soft repulsive part of the pair interaction potential. Along 
the MCT transition line in the Lennard- Jones liquid the temperature dependent packing fraction (/>(p, T) , resulting 
from the effective excluded volume diameter a = a°^{T) |6^, is roughly constant (f>{p,T) « 0.52 and approximately 
equal to its hard sphere value [Q. As the soft repulsion of the Lennard-Jones system leads to a"^ oc T^^/^^, the 
critical density smoothly decreases with temperature pof . Note that this observation indicates that the nonergodicity 
transitions resulting from the solution of Eq. |^ for the Lennard-Jones system are dominated by the excluded volume 
effect, i.e. the primary peak of the structure factor Sq, as is also the case for the hard sphere system. Because the 
idealized MCT has been developed from approximations aimed at describing the connected physical mechanism, called 
cage- or back-flow effect, the quality of the mode coupling approximation is expected to be unaffected by temperature 
changes as long as these excluded volume effects dominate in Eqs. ^ and |^. 

The nonergodicity transition lines of the HCAY system on the other hand are additionally affected by the low 
wavevector fluctuations in the fluid structure factor, Sq for q 0, and by the increase in Sq at large wavevectors 
arising from the short-range nature of the attraction. The first aspect, which also occurs for longer range attractions, 
leads to nonergodicty transitions tracking the spinodal curve (see appendix). 

The short-range nature of the attraction causes the stronger localization of the particles, i.e. the shorter localization 
length r^, upon decreasing the temperature. Also the strong increase in the elastic shear moduli along the nonergodicity 
transition line occurs only for sufficiently short-range attractions as shown in Fig. |^. Again, the ability of the MCT 
Eqs. |l| and ^ to describe such local interparticle correlations is not known. Note, however, that the HCAY results are 
independent of the numerical parameters chosen. Clearly, theories aimed at long wavelength phenomena at the gel 
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transition cannot incorporate these variations of the elastic modulus as described by the MCT because it arises from 
local potential energy considerations. 

The asymptotic model, defined by Eq. ^ (and b ^ oo), which highlights the effects of strong short-range attractions, 
captures all aspects of the low density MCT nonergodicity transitions qualitatively and even semi-quantitatively. 
Furthermore, it clearly demonstrates that the gel transitions are not driven by long-range structural correlations. 
It can be expected that such an asymptotic model can be found for other theories of liquid structure with strong 
short-range attractive potentials, but the detailed predictions presented here rest upon the use of the MSA for the 
HCAY system. 

Based on this suggested interpretation of the MCT nonergodicity transitions, several features of the computed HCAY 
density-temperature diagrams agree qualitatively with experimental observations made on colloid-polymer mixtures 
and sterically stabilized suspensions P- p^p^|jT7| ] . First, the gel transitions appear to lie at lower temperatures than, 
but otherwise track, the freezing line when present. Second, for short-range attractions the gel transition can shift to 
comparable or higher temperatures than those required for gas-liquid phase separation. Third, the gel line does not 
show such a strong density dependence as the static percolation transition. 

We emphasize that this suggested interpretation of the colloidal gel transitions is based on an extension of the 
idealized MCT of the glass transition beyond the range its approximations were aimed at. Our speculation, however, 
can be decisively tested by dynamic light scattering experiments. Nonergodicity transitions within the MCT exhibit 
universal dynamical properties p7|j2l| , |6l| , |6^ , which for example led to the identification of the colloidal hard sphere 
glass transition by van Megen and coworkers |^^|3^. As more complicated bifurcation scenarios, Ai with I > 2 
p7|,p9 40 6^, can be expected, the dynamics at the gel transitions should be very nonexponential and anomalous. 
Moreover, the short-range attractions lead to couplings among more wavevector-modes, as can be seen from the 
asymptotic model defined by Eq. |^, resulting in an MCT exponent parameter (see ref. for a definition and 

details on its calculation) A = 0.89 for b oo, considerably larger than values found for systems not characterized 
by short-range attractions (see e.g., pO|j4^ , |6^ ) . 

The proposed connection between the MCT nonergodicity transitions in the HCAY system and the non-equilibrium 
transitions in colloidal suspensions is further supported by the following observations. For moderate-range attractions, 
like the b=7.5 diagram in Fig. ^ the coexistence region can be tentatively divided into three regions. For somewhat 
lower temperatures than the critical temperature, gas-liquid phase separation occurs, provided a thermodynamically 
stable liquid phase exists. For temperatures (just) below the triple point temperature, gas-crystal phase separation 
takes place. Decreasing still more, gas-glass coexistence may be expected if — as argued from computer sim- 
ulations ]70[ | — the way to crystallization proceeds via the initial formation of a liquid droplet, whose density lies 
above the nonergodicity transition line. As the glass states for this system are rather close to the ideal hard sphere 
glass state, we expect signatures of this well studied transition to be observed [ p8| , p9pl| -|3^ . We speculate that the 
vanishing of the homogeneously nucleated crystallites in the colloid-polymer systems upon addition of sufficient large 
molecular weight polymer, observed in signals the presence of a nonergodicity transition as found in the colloidal 
hard sphere system |3jj7l|,0. This suggestion can be tested by st udying th e dynamic density fluctuations close to 
the transition as has been demonstrated in the hard sphere system [ P8| , p9y3l| -|34|] . 

For suspensions with short-range attractions, like the b=20- or b=30-curve in Fig. H, it seems possible that the 
long-range density fluctuations, likely induced by the hidden critical point, become arrested when the denser domains 
of the system cross the MCT nonergodicity transition line. Nonergodic gel states characterized by large small- 
wavevector form factors fq, rather short localization lengths, and finite, rather large elastic moduli can be expected. 
We suggest that these nonergodicity transitions cause the gel transitions observed in the colloid-polymer mixtures 
and sterically stabilized suspensions, and anticipate that they may also play a role in other colloidal systems, such as 
emulsions |73|, emulsion-polymer mixtures | |7^ , and suspensions of globular proteins p^-[25|], in which short-range 
attractions also dominate. We caution again that a proper extension of the MCT used here to include a full description 
of the critical dynamics close to critical points has yet to be formulated. Experimental tests of the dynamics close to 
the gel transitions would be required to test our suggestion. 
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APPENDIX: 



In this appendix the nonergodicity transitions caused by the increase in the q hmit of the structure factor close 
to the spinodal hnes are discussed for the HCAY system. Figure shows the spinodal hnes and the gel transitions for 
the interaction parameters considered in Sections IV and Also shown are nonergodicty transition lines occurring 
only close to the spinodal lines. As seen, there is at least one crossing of the two types of nonergodicity transitions for 
each attraction range b, where in all cases the gel transitions discussed in the main text provide the larger, physical 
nonergodicity parameters. The additional transition lines presented here in the appendix have two peculiarities which 
cause us to doubt the validity of the present MCT for their description. First, they are directly caused by the small 
wavevector structure in Sq; thus, a proper MCT description should include also the very likely present critical dynamics 
p3| . Second, at these transition lines only the collective density fluctuations for exceedingly small wavevectors or on 
large length scales are arrested. The single particle dynamics remain fluid-like, i. e. /| = from Eq. 2. Again, this 
suggests that long-range collective fluctuations are of crucial importance at these transitions and the simple MCT 
approach used here is likely insufficient. 
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FIG. 1. AHS phase diagram in terms of the reduced temperature r and the particle volume fraction 4>. The labeled lines are 
the bifurcation lines of Eq. |l| discussed in the text. The shaded region encloses nonergodic density fluctuations. Liquid-glass 
transitions occur when crossing the Bl and B2 lines into the shaded nonergodic region. Glass-glass transitions occur along 
the piece of the B2 line inside the nonergodic region. The physically relevant endpoint singularity is indicated by a circle. The 
data for the fluid-solid freezing and melting lines (•) are taken from MWDA calculations |55|. The arrows show the (PY) hard 
sphere freezing (0f) and melting (0m) volume fractions 

FIG. 2. HCAY phase diagrams in terms of the reduced temperature and particle volume fraction (j) for varying values of 
the screening parameter b, as labeled. MSA gas-liquid spinodals are shown with critical points denoted by (•), together with 
the corresponding MCT glass transition lines as labeled. The chain curves correspond to the asymptotic prediction in Eq. ^ 
with Fcffe oo) = 3.02. The open circles (o) denote the locations in the b=20 diagram of the suspensions for which the results 
in Fig. H were calculated. 

FIG. 3. HCAY critical coherent nonergodicity parameters fq for b=30 as functions of normalized wavevector qa and Yukawa 
prefactor K. Note that the volume fraction varies according to the b=30 critical glass transition line in Fig. ^ and that the 
K-values jump from K = 14 to 20 and 50. Rear curve shown in bold is the density and temperature independent asymptotic 
prediction resulting from Eq. ^ and b oo. 

FIG. 4. Localization length (root-mean-square displacement) in the glass along the critical glass transition lines in Fig. ^ 
as a function of the reduced temperature for various values of the Yukawa screening parameter b, as labeled. 

FIG. 5. Elastic shear modulus (in units of kT/a^) in the glass along the critical glass transition lines in Fig. ^ as a function 
of the reduced temperature K^^ for various values of the Yukawa screening parameter b, as labeled. 

FIG. 6. Normalized intermediate scattering function at a fixed wavevector qcr=Q.2 as a function of Yukawa prefactor and 
dimensionless time, as labeled. The corresponding locations in the b=20 HCAY diagram are shown in Fig. ^ as open circles. 

FIG. 7. HCAY phase diagrams showing, in addition to the MSA gas-liquid spinodal curves and the MCT glass transition 
lines as in Fig. |2|, the MCT transition lines (short-dashed lines) which are characterized by a small-g structural arrest in fq 
and fq = 0. The inset shows an enlargement of the line crossing in the vicinity of the b=20 critical point. 
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